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■ In this work, within the QCD factorization approach, we investigate the branching frac- 
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Os| ■ tions and CP asymmetries of decays B -> £q(1430)p and B -> Kfi (1430)0 under two dif- 
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ferent scenarios both in the standard model and the family nonuniversal Z' model. We find 
that the annihilation terms paly crucial roles in these decays and lead to the main uncer- 
tainties. For decays B — > Kq (1430)p°(ffl), the new Z' boson could change branching 



! fractions remarkably. However, for other decays, its contribution might be clouded by large 

uncertainties from annihilations. Unfortunately, neither the standard model nor Z' model 
can reproduce all experimental data under one certain scenario. We also noted that the CP 
asymmetries of B — > (1430)p°(to) could be used to identify the ^(1430) meson and 
| search for the new physics contribution. 

o 

O I. INTRODUCTION 



The study of B meson rare decays is a crucial tool in testing the fundamental interactions among 
elementary particles, exploring the origin of CP violation, and searching for possible new physics 
(NP) beyond the standard model (SM). Theoretically and experimentally, such kind of research has 
been conducted in great detail, especially in the weak interactions of B meson. In particular, the 
processes induced by flavor-changing neutral-current (FCNC) only occur at the loop level in SM, 
and are therefore a very sensitive probe of NP beyond SM . Already, FCNC processes have been 
explored mainly in the B q — B q mixing and the semi-leptonic weak decays, which permit a clean 
theoretical description. So far, the charmless hadronic B meson decays induced by FCNC have 
also been studied extensively, such as B — > Kn^K^'fy and K^rj^ decays. In the past few years, 
the new physics effect in these decays have also been studied widely, such as in supersymmetry 
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model, two-Higgs doublet model, Z ! model, the forth generation model, extra dimension models, 
and so on (see review in [1] and references therein). 

In order to search for effect of NP in the nonleptonic B decays, most theoretical studies are 
focused on B — > PP, PV or VV in the past few years. But, the studies of decay modes involving 
a scalar meson are relatively few, because the underlying structure of the scalar mesons is not 
well established in theoretical side. To describe the component of the scalar mesons, there are 
usually two possible scenarios (SI and S2) according to the QCD sum rule method [2]: (i) In 
SI, we treat scalars above 1 GeV as the first excited states, while the scalars under 1 GeV are 
regarded as the low lying states; (ii) In S2, the scalars above 1 GeV are viewed as the ground 
states, and light scalars are four-quark bound states or hybrid states. Under these two scenarios, 
many special decays have been examined within the QCD factorization (QCDF) approach [3, 4] or 
the perturbative QCD approach (pQCD) [5-14]. However, because of large uncertainties in SM, 
the NP effects in these decays are rarely studied. 

Very recently, BaBar collaboration reported their first branching fraction measurements for the 
decays B ->■ K*{U30)p that are induced by FCNC [15]: 

Br(B° '(1430) p°) = (27 ±4±2 ±3) x 10~ 6 ; (1) 

Br(B° ^K* + (U30)p-) = (28 ± 10±5±3) x 10~ 6 . (2) 

The above results are inconsistent with the pQCD predictions [12] in most cases. Moreover, these 
results are somewhat much lower than the QCDF predictions [4] but are consistent with QCDF 
within rather large uncertainties. For B — > Kq (1430)0, BaBar collaborator also updated their 
results [16, 17] in ref. [18]: 

Br(B°^K£°( 1430)0) = (4.3 ±0.6 ±0.4) x 10~ 6 ; (3) 
5^^^(1430)0) = (7.0±1.3±0.9) x 10~ 6 . (4) 

Both QCDF and pQCD calculation of above modes have also been presented in Refs.[4, 8], and 
the predicted central values of 5° — > 0^q°(143O) deviate from the experimental data, though they 
can be also accommodated within very large theoretical errors. In the following, Kq * (1430) is 
denoted as in some places for convenience. 

The predictions of SM cannot agree the data convincingly, which gives us possible hints on 
physics beyond SM. It is our purpose of this work to show that a new physics effect of similar size 
can be obtained from some models with an extra spin-1 Z' bosons, which are known to naturally 
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exist in some well-motivated extensions of the SM [19]. Interesting phenomena arise when the 
Z ! couplings to physical fermion eigenstates are nondiagonal, which could be realized in the E 6 
models [20], string models [21] and some grand unified theories [22]. For example, in the super 
string model advocated by Chaudhuri et.al.[2l], it is possible to have family nonuniversal Z' cou- 
plings, because of different constructions of the different families. It also should be note that in 
such a model, called the family nonuniversal Z' model, the nonuniversal couplings could lead to 
FCNCs at the tree level as well as introduce new weak phases [23], which could explain the CP 
asymmetries in the current high energy experiments. In fact, the effects of Z' models have been 
studied extensively in the low energy flavor physics phenomena, such as neutral mesons mixing, 
B meson decays, single top production and lepton decays [23-31]. 

In this current work, we shall adopt the QCD factorization approach [32] to evaluate the relevant 
hadronic matrix elements of B decays, since it is a systematic framework to calculate these matrix 
elements from QCD theory, and holds in the heavy quark limit m b — > °o and the heavy quark 
symmetry. In such calculations, one requires the additional knowledge about form factors of B 
meson to the scalar or the vector transitions. This problem, being a part of the nonperturbative 
sector of QCD, lacks a precise solution. To the best of our knowledge, a number of different 
approaches had been used to calculate the form factors of B — > S decays, such as QCD sum rule 
[33, 34], light-cone QCD sum rule [35, 36], perturbative QCD approach [37] and covariant light 
front quark model (cLFQM) [38]. Among them, the form factors of the cLFQM are first calculated 
in the spacelike region and their momentum dependence is fitted to a 3-parameter form. This 
parameterization is then analytically continued to the timelike region to determine the physical 
form factors at q 2 > 0. Moreover, for these form factors both the heavy quark limit and heavy 
quark symmetry are satisfied. For that reason, we will use the results of cLFQM [38] in the 
following calculations. 

For comparison, B — >■ K^p and KqQ decays in SM should be reinvestigated in Section.II. In 
Section. Ill, we will review the family nonuniversal Z' model briefly and show the effect of Z' 
to decay modes we are considering. In Section.IV, we will present our numerical results and 
discussions in great detail. At last, we will summrize this work in Section. V . 
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II. REVISITING B -> pK* (1430) AND B -> <j)K* (1430) DECAYS WITHIN THE QCDF 

FRAMEWORK 

To proceed, we discuss the decay constants of the scalar meson. Unlike pseudoscalar meson, 
each scalar meson has two decay constants, the vector decay constant f$ and the scale-dependent 
scalar decay constant f$ namely, which are defined as: 

(S(p)\q 2 Ynqi\0) =fsPn, (S(p)\q 2 qi\0)=msfs, (5) 

and they are related by the equation of motion: 

m 2 (n)-mi(n) - 

JS = Js, (6) 

ms 

where m 2 and m\ are the running current quark masses. Therefore, the vector decay constant is 
much smaller than the scalar one. As for the vector meson, the two kinds of decay constants are 
also given by [39] 

(V( P )\q 2 7nqi\0) =fvm v e;, (V(p,E*)\qa^q'\0) = /^e* - Pv e*). (7) 

The twist-2 and twist-3 light-cone distribution amplitudes (LCDAs) of scalar mesons, <j>s(x), 
and <j>°(x) respect the normalization conditions: 

/»' "* w = r w = /.' = w 6 ' <8) 

and (f)J(x) = ^-^^s( x )- The twist-2 LCDA can be expanded in the Gegenbauer polynomials: 

5 (x,ju) = -^/s(/i)6x(l -x) £ 5 m (ju)d /2 (2x- 1). (9) 
v6 m=1 

The decay constants and the Gegenbauer moments of the twist-2 wave function in two different 
scenarios have been studied explicitly in Refs. [3] using the QCD sum rule approach. As for 
the explicit form of the Gegenbauer moments for the twist-3 wave functions, there exist some 
uncertainties theoretically [40], thus we choice the asymptotic form for simplicity: 

*i=^//» *J = ^/5(l-2x). (10) 

For the vector mesons, the normalization for the twist-2 function <Ity and the twist-3 function 
4> v is given by 

! dx<f> v (x) = f v , [ dx® v (x)=0, (11) 
Jo Jo 
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where the definitions for <5 v (x) can be found in [32]. The general expressions of these LCDAs 
read 



<*V(x,ju) =6x(l-x)/ v 



l + £a:( M )C„ 3/2 (2*-i; 

n=l 



(12) 



and 



2^-l + £a„ v _ L ( A i) J P„+i(2^-l) 



(13) 



where P n {x) are the Legendre polynomials. 

In the calculation, the most important nonperturbative parameters are form factors of B — > 5, V 
transitions, which are defined by [41]: 



(V(p')\V,\B(p)) = --^—e^ £ * v P a g Pv BV ( q 2 ), 

niB + my n 



(V(p')\A^\B(p)) = i\(m B + m v )e;A B l v (q 2 ) ^-P M Af (q 2 ) 



m B + m v 

-2m,^,[AfV)-An<7 2 )]}, 



(S(p')\A^\B(p)) = -i 



2 2 

m B~ m S „ \ Z7S5/ 2n , m B~ m S„ j?BS(Jl\ 



^ Kl (^ ) + 



(14) 



withfy = (p + p')n,q^ = (p-p')n- 

To calculate the amplitudes, we start from the effective Hamiltonian responsible for & — >■ 5 
transitions, which is given by [42] 



G F 



10 



Vutv* (Ci o? + c 2 02) + v cfe v* (d of + c 2 o|) - v,*v* ( £ qo, 

i=3 



V2 

+ CqyOly + CggOsg) 



+h.c. 



(15) 



In the above equation, V q bV* s (q = u,c,t) represent for products of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements, Q are the responding Wilson coefficients, and 0[ are the rele- 
vant four-quark operators whose explicit forms could be found, for example, in Refs. [42]. 

We now turn to study the short-distance contributions within the QCDF approach, where the 
contribution of the nonperturbative sector is dominated by the form factors and the nonfactorizable 
impact in the hadronic matrix elements is controlled by hard gluon exchange. The hadronic matrix 
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elements of the decay can be written as 

(MiM 2 \O t \B) = J^lf*" 1 J* thTfj(x)* Mi (x) 



j 

+ [ l dt[ l dx ['dyT/'i^x^B^M^Mziy), (16) 
Jo Jo Jo 

where T/j and T/ 1 denote short-distance interactions and can be calculated perturbatively. <t>x(x) 
are the universal nonperturbative light-cone distribution amplitudes. Using the weak effective 
Hamiltonian given by Eq.(15), we then obtain the decay amplitudes as: 

V2 p=U;C [\ ll j ^ 

x 2/^^(fnJ)mBPc-/B/^(W + fc3 + ^)^J, (IV) 

A(B° -^K* Q °(j)) = & £ ^P S) \ ( a 3+ a 4+ a 5-4( a 6-l a s)-l( a T + a 9+ a lo)) 
V Z p=u,c [\ Z Z J 

x 2f^ Ko (ml)m B p c -f B fyf K *(b 3 -^b 3) Ew) K ^y (18) 



A(B-^T p-) = £L ^ { " (< + r X « ~ \<) ~ ^io) 



P K o 



x 2f K *AQ p (ml*)m B p c -f B f p f KS (b 2 8!? + b 3 +b 3 , EW ) pK » j>, (19) 



x 2f K *A Q IJ (m K *)m B p c + 



3 

+ -(a 9 + a 7 ) 



IfpF™ (m 2 p )m B p c 



Hp 

- f B fpf K *(b28P + b 3 + b 3m ) pK ^, (20) 

A(B°^^-p+) = i% £ - («i^ + a£ + H?«g + «f + r?af) 

V2.p= u ,c { v y P^o 



< 2f K *AQ p (m 2 K *)m B p c - f B f p f K * (b 3 - ^b 3 , E w) pK * \ , (21) 



1 
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F*0 o 



p=u,c 
i B P/^2 



-a p 4 -4°(a p -^a p ) + ^a p 



1 



1 



X 2fKZ A ( m K*) m BPc + 



a 2 8 p + -(a 9 + a 7 ) 



2f p Fi K °(m 2 )m B p c 



- fBfpf K * {~b 3 + ^3,ew) pK * } , 



(22) 



p=u,c 



•p 



(32^ + 2(a3+a 5 ) + -(a9 + a 7 ) 



2fcoFi K ° {m 2 a )m B Pc 



K* Q G) 



- (a x 8 P +a p A + rj° (aj + of) + af ) ^ 2/^A^ ffl (m|*)m g p c 

- fBf(ofK*(b 2 8P + &3+^3,Ew) ^ f, 



(23) 



A(5 ->^ co) 



4 I 

p=u,c 



a 2 5jf + 2(a 3 +a 5 ) + -(a 9 + a 7 ) 



- (< + rJ°(a£--a£)--af ) 2/ z *Ao P (m|* )m g p 6 
z z / raz* 



2/coFj °(m Q) )m B p c 



~ fBfcofK* (h ~ 2^3,Ew) 



(24) 



= 



2m| 



(25) 



where the ratios r!^ and r£ are defined as 

* m fo (/i) / v ' '* VA ~ ; m z ,( J u)(m 2 (ju)-mi( J u))' 

The order of the arguments of the af (MjM 2 ) and bi{M\M 2 ) coefficients is dictated by the subscript 
M\M 2 , where Mi shares the same spectator quark with the B meson and M 2 is the emitted meson. 
For the annihilation part, M\ is referred to the one containing an anti-quark from the weak vertex, 
and M 2 contains a quark from the weak vertex. Combining the short-distance nonfactorizable 
corrections, the effective Wilson coefficients af have the expressions 



a p {M x M 2 ) 



c i + %lW 2 ) + C ' ±lCFa * 



An 1 

V l (M 2 ) + —H i (M l M 2 ) +P?{M 2 ), (26) 



N c J N c 4n 

where Vi(M 2 ) account for vertex corrections, Hi{MiM 2 ) for hard spectator interactions and P,(M 2 ) 
for penguin contractions. The coefficients b\ and bi^w stand for the contribution of annihilation 
diagrams. 
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In QCDF approach, the end-point singularities appear in calculating the twist-3 spectator and 
annihilation amplitudes. Since the treatment of endpoint divergences is model dependent, sub- 
leading power corrections generally can be studied only in a phenomenological way. As the most 
popular way, the end-point divergent integrals are treated as signs of infrared sensitive contribu- 
tions and parameterized by [32]: 

/'^Ux^d+p^ln^, ,27) 
Jo y A h 

with the unknown real parameters pA and 0a- More discussion about them will be in Section.IV. 

HI. THE FAMILY NONUNIVERSAL Z' MODEL 

In this section, we will review the main part of the family nonuniversal Z' model briefly. In 
the current work, for simplicity, we only focus on the models in which the interactions between 
the Z' boson and fermions are flavor nonuniversal for left-handed couplings and flavor diagonal 
for right-handed cases. Of course, the analysis can be straightly extended to general cases in 
which the right-handed couplings are also nonuniversal across generations. The basic formulas 
of the Z' model with family nonuniversal and/or nondiagonal couplings have been presented in 
Refs.[19, 23], to which we refer readers for detail. Here, we just review the ingredients needed in 
this work. 

In the gauge basis, the neutral current Lagrangian induced by the Z' boson can be written as 

= -82^ , (28) 

where g2 is the gauge coupling associated with the additional U(\)' group at the Mw scale. Ne- 
glecting the renormalization group (RG) running effect between Mw and M z i and the mixing 
between Z' and Z boson of SM, we present the chiral current as 

^ =YMy* [MjPl + ^jPr] Wj , (29) 

hi 

where the sum extends over the flavors of fermions, the chirality projection operators are P^r = 
(1 T Ys) /2, the superscript / stands for the weak interaction eigenstates, and £y L {£y R ) denote the 
left-handed (right-handed) chiral couplings. e VL and £y R are required to be hermitian so as to 
arrive a real Lagrangian. Accordingly, the mass eigenstates of the chiral fields can be defined 
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by Vl,r = VvlrWlr> an ^ tne usua l CKM matrix is given by Vckm = Vu L Vj L - Then, the chiral Z' 
coupling matrices in the physical basis of up-type and down-type quarks are, respectively, 

B X U = V ux e ux vl , fl} = V dx E dx vl (X = L,R). (30) 

If the £ matrices are not proportional to the identity, the B matrices will have non-zero off-diagonal 
elements, which induce FCNC interactions at the tree level directly. In this work, we assume that 
the right-handed couplings are diagonal for simplicity. Thereby, the effective Hamiltonian of the 
b — > sqq(q = u,d) transitions mediated by the Z' is 

■*df = ^| (f^) B^(b S ) v ^(B L qq (qq)v-A + B R qq (qq) v+A ) +h.c. , (31) 

where g\ = e/ (sinOiy cos 6w) and M z t the mass of the new gauge boson. We note the above oper- 
ators of the forms (bs)y-A(qq)v-A and (bs)v-A(<jq)v+A already exist in SM, so that we represent 
the Z' effect as a modification to the Wilson coefficients of the corresponding operators. Hence, 
we rewrite the eq.(31) as 

Jtlg = (aC 3 0^ + AC 7 0\ q) + ACgO^) +h.c, (32) 

where the additional contributions to the SM Wilson coefficients at the M w scale in terms of Z' 
parameters are given by 

Thus we can have a Z' contribution to the QCD penguins AC 3 ( 5 ) as well as the EW penguins 
AC 9 ( 7 ), in the light of the results found by Buras et al. [42]. In order to show that the new physics 
is primarily manifest in the EW penguins, we assume B^^ — which have been used 

widely [25, 27, 28, 30]. As a result, the Z' contributions to the Wilson coefficients at the weak 
scale are 

AC 3(5) =0, (35) 



AC 9{1) =4 l ^^e-^, (36) 
v tb Vts 
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where 



lx = ( giM z \ 2 



\g\M z >) 
fa = Arg[5^] . 




(X = L,R), 



(37) 



(38) 



Because of the hermiticity of the effective Hamiltonian, the diagonal elements of the effective 
coupling matrix must be real. However, the off-diagonal elements, such as B^ b , generally may 



of universality for the first two families, as required by K and ji decay constraints [23]. Since 
the major objective of our work is searching for new physics signal, rather than producing acute 
numerical results, we also assume B qq ~ B qq , because we expect that \B qq \ and \B qq \ should have 
the same order of magnitude. 

It should be emphasized that the other SM Wilson coefficients may also receive contribu- 
tions from the Z' boson through renormalization group (RG) evolution. With our assumption 
that no significant RG running effect between M' z and M\y scales, the RG evolution of the 
modified Wilson coefficients is exactly the same as the ones in SM [42]. The numerical re- 
sults of Wilson coefficients in the naive dimensional regularization (NDR) scheme at the scale 
/l = 2.1GeV (Hh = lGeV) are listed in Table I for convenience. 

In summary, we list here our simplifications to a general Z' model: we assume (i) no right- 
handed flavor-changing couplings (Bfj = for i ^ j), (ii) no significant RG running effect between 
M z i and M\y scales, (iii) negligible Z' effect on the QCD penguin (AC3.5 = 0) so that the new 
physics is manifestly isospin violating, (iv) \B qq \ and \B qq \ are same so as to reduce the number of 
parameters. With these simplifications, we have only two parameters left in the model. So, this 
approach provides a minimal way to introduce the Z' effect in the concerned decay modes. Of 
course, more general Z' models are possible. 

Now, the only task left is to constraint the parameters within the existing experimental data. 
Generally, g2/gi ~ 1 is expected, if both the U(l) gauge groups have the same origin from some 
grand unified theories. We also hope Mz/M z t ~ 0.1 so that TeV scale neutral Z' boson could be 
detected at LHC. Theoretically, one can fit the left three parameters \B^ b \, \Bj d \ and new weak 
phase fa with the accurate data from B factories and other experiments such as Tavatron and LHC. 
For example, B^ b and fa could be extracted from B s -B s mixing as well as B — > K^£ + £~ decays. 
To resolve the mass difference between B s and B s , \B^ b \ ~ IV^V^I is required [25, 28, 43]. In 
Refs.[28], the authors got the fa is about -80° by fitting data of B s -B s mixing and B -> K^l + T 
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TABLE I: The Wilson coefficients C, within SM and with the contribution from Z' boson included in NDR 
scheme at the scale jx = 2.1 GeV and \ih = 1.0 GeV. 



wuson 
coefficients 


jU = 2.1 GeV 


Hh = 1.0 GeV 


i^SM 


ACf 


i^SM 


ACf 


Ci 


1.135 





1.224 





C2 


-0.283 





-0.429 





C3 


0.021 


0.09^-0.02^ 


0.034 


0A5c, LL - 0.04^ 


c 4 


—0.049 


r\ in!: LL i a a i £ LR 
— 0.20g +O.OI5 


—0.072 


f\ 1 1 £ LL i n no £ L/? 

— 0.31§ + 0.03q 




0.010 


0.03g"' + 0.02^" < 


0.010 


r\ r\n f I T , r\ r\r* f f R 

0.02^ + 0.025^ 




—0.06 


a j; 7,7, , r\ ao e 7,7? 

— O.265 + 0.035 


—0.104 


— 0.445 +0.075 


/,-/ 


-0.018 


5.3^ LL -461^ LS 


-0.023 


6.3^-457^ 


Cs/ oc em 


0.081 


2.43^ LL -286^ LS 


0.134 


4.8^-497^ 


C9/ a em 


-1.266 


_ 594 £LL + 6A £LR 


-1.366 


-643^ LL + 7.8^ LR 


C10/ &em 


0.321 


m$ LL -\.o$ LR 


0.483 


257^ LL -1.9^ LS 


Cjy 


-0.345 




-0.395 






-0.161 




-0.181 





decays. Subsequently, with B L sb and 0l arrived and experimental data of B — > KK,Kn,Kp and 
k(*)$, Bq q and B^ q could be extracted analogously. Specifically, the CP asymmetries in B — > 

j Li R L R 

K§,Kk can be resolved if \B^ b B S s | ~ | V^^V^* | , which indicates \B qq | ~ 1. However, we have one 
remark here. In dealing with the nonleptonic B decays, because different groups used different 
factorization approach, the fitted results are different, but all results have same order. Noted that 
the detailed constraint of these parameters is beyond the scope of current work and can be found in 
many references [27, 28]. Summing up above analysis, we thereby assume that t, = ^ LL = ^ LR £ 
(1(T 3 , 1CT 2 ) and <fc e (-60°, -90°) so as to prob the new physics effect for maximum range. 

IV. NUMERICAL RESULTS AND DISCUSSION 

In this section, to begin with, we will give the parameters used in this work. Since it is not clear 
whether the scalar meson Kq (1430) belongs to the first orbital excited state (SI) or the low lying 
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resonance (S2), we will calculate the processes under both scenarios. In the calculation, the decay 
constants and Gegenbauer moments obtained within the QCD sum rules method under different 
scenarios are presented as follows [3]: 

SI : /^(l.OGeV) = -300MeV;/^(2.1GeV) = -370MeV;Bi(1.0GeV) = 0.58; 

fii(2.1GeV) = 0.39;5 3 (1.0GeV) = -1.20;fl3(2.1GeV) = -0.70; (39) 
S2 : /^(l.OGeV) = 445MeV;/ z *(2.1GeV) = 550MeV;fli(1.0GeV) = -0.57; 

Bi(2.1GeV) = -0.39;fl 3 (1.0GeV) = -0.42;fl 3 (2.1GeV) = -0.25. (40) 

In QCD sum rules method, the major parameter is the Borel window, which takes large uncertainty 
to the parameters listed above. In Ref.[3], the authors had discussed the errors caused by them in 
great detail and found that #13 will take 30% changes. As a result, we will not discuss this part 
any more in the current work. 

For the vector mesons, the longitudinal and transverse decay constants are list as: 

f p = 216 MeV, f a = 187 MeV, / = 215 MeV, 

fp 1 = 165 MeV, /i = 151 MeV, f£ = 186 MeV, (41) 

where the values are taken from [44]. In the LCADs of vectors, the Gegenbauer moments a„ and 
a n± nave b een studied within the QCD sum rule method. Here, we will employ the most recent 
updated values [45] 

a£' ffl = 0.15, ajf = 0.14, of =0.18, a| ± = 0.14, (42) 
and a\ = 0, af ± = 0. 

As stated earlier, various form factors for B — >■ 5, V transitions have been evaluated in cLFQM 
[38]. In this model form factors are first calculated in the spacelike region and their momentum 
dependence is fitted to a 3-parameter form 



F(0) 

\-a{q 2 /m 2 B )+b{q 1 /m 2 B f 



The parameters a, b and F(0) relevant for our purposes are summarized in Table.II. 

In Refs.[3, 4, 14], it was found that in decay modes with scalars the main theoretical uncertain- 
ties are due to the weak annihilations, especially for the penguin dominated ones. In B — > PP,PV 
decays, the annihilation amplitudes are helicity suppressed because the helicity of one of final 
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states cannot match with that of its quarks. However, this helicity suppression can be alleviated 
in the decay modes with scalar because of nonvanishing orbital angular momentum. Thus, weak 
annihilation contribution to B — > SP(V) is much larger than the B — > PP(V) case. However, as 
stated before, the end-point singularity appears in calculating the annihilation contribution, and 
then two free parameters, p A and <j> A , are introduced phenomenally. In Ref. [4], it is found that 
the behavior of SV is similar to the longitudinal part of VV. Fortunately, with experimental data, 
it presents the moderate value of nonuniversal annihilation phase (j> A = —40° for B — y VV decay 
modes [32]. Therefor, for B ->■ SV, we conservatively take (j) A = (-40±20)° with p A = 0.6 ±0.2, 
which also assures that the hadronic uncertainties are considerably reduced. Furthermore, the 
endpoint divergence Xh in the hard spectator contributions can also be parameterized in the same 
manner. 

Within above parameters and formulas, we calculate the branching fractions of these decays 
in SM and the family nonuniversal Z' model under two different scenarios. Together with partial 
experimental results, the results under are exhibited in Table.III, respectively. For the center values, 
we adopt t, = 0.005 and (j)[ h = —80°. For all theoretical predictions, the first errors arise from 
the power corrections of weak annihilation and hard spectator interactions characterized by the 
parameters X A h- To obtain the second errors of the Z' model results, we scan randomly the points 
in their own possible parameter spaces. 

Comparing our predictions of SM with those in Ref. [6] (considering the typos), there are few 
differences. Some reasons are list as follows: (1) In the Ref. [6], for the parameterizations of 
singularities, the center values correspond to p AH = and § a ,h = 0, while we set p AH = 0.6 and 
<j> Aj H = —40°; (2) The difference of Wilson coefficients, caused by the top quark mass and other 
part parameters, will change the results slightly; (3) In this work, the different form factors of 
B — > ^(1430) are used under different scenarios, but they adopted same values in the ref. [4]. 

TABLE II: Form factors of B — > p,^ ( ,(1430) transitions obtained in the covariant light-front model [38]. 



F F(0) F(q 2 m J a b 


F F(0) F(q 2 max ) a b 


V B P 0.27 0.79 1.84 1.28 
Af 0.22 0.53 0.95 0.21 
Ff^tSl] 0.21 0.52 1.59 0.91 
*f*"[S2] 0.26 0.70 1.52 0.64 


Aq P 0.28 0.76 1.73 1.20 
A B 2 P 0.20 0.57 1.65 1.05 
F B *°*[S1] 0.21 0.30 0.59 0.09 
F BX °*[S2] 0.26 0.33 0.44 0.05 
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TABLE III: Branching fractions (in units of 10 6 ) under the different Scenarios 



Decay Mode 


SI 

SM SM+Z' 


S2 

SM SM+Z' 


Expt 


B ^K*-<p 
5° -> K*°(j) 
B ^K*°p- 
B ^K*~p° 

r 

B° -»• ^°p° 

B ^K* co 
B° -»• ^°0) 


9 4+5.1 o 0+5.6+1.8 
Z - H -1.8 -'•°-2.5-1.2 

9 9+ 4 - 9 4 7+5.2+4.4 
z,.z._ 17 ^- '_2.6-2.2 

1 1 7+ 8 - 4 1 1 ^+7-7+3.6 
11 -'-4.0 11 - J -4.7-3.8 

7 7+4.5 107+5.9+27.4 
—2.3 —3.4—10.0 

4 6 +L7 3 o+i - 4 + 6 - 4 

H - D -0.7 -'•*-0.8-1.8 

10 7+ 8 5 l4 4+ 9 - 2 + 8 - 7 

lU./_ 37 l^-^_5.i_4.i 

^ 6+ 3 - 1 7 Q+3.7+10.1 
J-D_14 '-0-1.7-3.8 

J -*-0.6 H - u -().7-0.6 


22 6 +19J 16 5+18-5+5.0 

zz - D -8.6 1D - J -10.3-7.9 

77 4+19.4 71 7 + 18.9+3.0 
zz - H -8.4 zi - z -8.4-4.6 

45 5+20.6 4, 4+19.7+4.4 

17.6+2° 15.9+*,+^ 

—4.4 —6.1—5.6 

94 <5+ 7 -° o+7.7+32.1 
ZH - J -3.8 - , -'- :, -4.8-8.5 

44 7+21.3 ^4 1+22.4+25.7 
^•'-10.5 JH - 1 - 12.6-8.6 

12.6™ 13.7 ^ 
10.6+1? 10.71^^ 


7.0±1.3±0.9 
4.3 ±0.6 ±0.4 

27 ±5.5 

28 ±12 



In Tables. Ill, for Kq (1430)0 channels, though the central values of the predicted under SI are 
smaller than the experimental data, they are accommodated with the large uncertainties. However, 
in S2, the theoretical results are much larger than the data, and cannot agree with data even with 
uncertainties. It should be noted that in this work we have not included the errors from uncertain- 
ties of meson distribution amplitudes (B\ and B3). Even with those uncertainties, the theoretical 
results are still larger than the upper limits of the data. These theoretical results also agree with 
the results from pQCD approach [8]. For 5° — y Kq°p° channels, contrary to £q (1430)0, the re- 
sult of S2 agree with data well and the prediction of SI is much smaller than the data. Since for 
B° — y Kq~p + there is large uncertainty in the experimental data, the theoretical results under both 

51 and S2 can accommodate the data with large uncertainties theoretically. That's to say, it is 
impossible to explain all data under one settled scenario simultaneously. 

When adding the contribution of the Z' gauge boson, as shown in the table, the Z' gauge bo- 
son changes the branching fractions under both two different scenarios. For B — y A^0 channels 
dominated by the weak annihilation, the Z' will enhance the branching fractions in SI, while in 

52 the branching ratios are decreased. The reason is that the weak annihilation is proportional 
to the decay constant f K *, which has different sign in different scenarios. For B~ ^K*°p- and 
5° — y Kq~ p + , as the scalar particle is the emitted particle, the whole amplitudes are proportional 
to the decay constant f K *, thus the new physics contribution have same behavior in different see- 
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narios. For channels with p° or CO, the spectator quarks enters not only the scalars but also the 
vectors, the amplitudes become more complicate, and we cannot describe the relation between 
new physics and branching fractions apparently. 

Compared to the experimental data, the Z' boson could change the branching fractions remark- 
ably and alleviate the disparities. However, we cannot achieve a definite conclusion yet whether 
Kq belongs to the ground states or the first orbital excited states. Moreover, for most modes except 
B ^K*°-p°(co), the new physics contribution might be clouded by the uncertainties taken by 
the weak annihilations. Thus, it is also very difficult to search for Z' effect in these decays. Specif- 
ically, for decays B — > Kq~p°(oj), Z' boson could enhance the branching fractions more than 2 
times, we hope these two channels could be measured in the LHC or Super-b factories in future so 
as to probe the Z' gauge boson. 

To test the isospin symmetry and prob new physics, we define two ratios: 

Br(B° -*K*°p°) +107 r 

* = Mg°^fp+r ft% -° : " |ExpJ; <44) 

t(B~) Br(B° ->• K*°<j>) -°- 51 1 FJ 

where the experimental results are also given and all uncertainties are added in quadrature. In the 
isospin limit, Ri = 1/2 and R 2 = 1 are expected to hold. Here, we list the theoretical results under 
different scenarios in different models: 

' n 97+0.15+0.31 <;,. 
u - DZ -0.09-0.07' ,5Z - 

. j 1.00+g-gJ, SI; r h I 0.76+Jj^, SI; 

R2[SM] = \ ; i? 2 [SM + Z'] = ^ To To I (47) 

n Q/i+u.oo 09 a 70+0.10+0.19 09 

In the above results, the theoretical uncertainties are reduced since they are ratios of branch frac- 
tions. We see that the symmetries are almost held in SM. However, the data shows that the isospin 
symmetries are violated, which means that the large weak annihilation may break the isospin sym- 
metry remarkably. When adding Z' contribution, except Rl under S2, the isospin symmetries are 
broken in an opposite direction. However, the family nonuniversal Z ! model cannot be ruled out 
due to large uncertainties in the experiments. 

Finally, we will discuss the CP asymmetries of these decays. For the charged mode B — > 
Kq~§, because \V u bV us \(?i 4 ) <C IV^V^I (A 2 ) and there is no tree contribution in the neutral mode 
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TABLE IV: The direct CP asymmetry (%) under the different scenarios 



Decay Mode 


SI 

SM SM+Z' 


S2 

SM SM+Z' 


B -y Kq°p ~ 
B -»• K*~p° 
B°^K*°p° 
B Q ^K*-p + 
B~ —y Kq~ (0 
B° -> f*°0) 


4+4 o+2+4 

3 1 Z 

q+26 9A +52+21 
^-38 ZH -42-12 

i+l o+l+l 
L -2 ^-1-0 

o+6 4 +2+5 
J -7 ^-3-3 

1 6+ 26 1 7+28+7 
10_ 39 l/_ 40 _5 


9+2 9+2+O 
— 1 -1-0 

1+3 £+3+10 
t — J — o 

ii+lO q+9+4 
u -13 v 112 

i +0 i +0+0 
L -0 1 -0-0 

i +4 4+3+5 
l -5 ^-4-6 

1Q+14 4+15+19 
ly -15 ^ 19 13 



B — y Kq ^), the direct CP asymmetries are almost zero in both SM and the Z' model. For B^K^p, 
although the CKM elements are suppressed, the tree operators with large Wilson coefficients ap- 
pear in the emission diagrams, so the amplitudes of tree and penguin may have comparable mag- 
nitudes. Thus, large CP asymmetries in these decays are expected, just like decays B — y Kit and 
B — y Kp . In Table.IV, we give the CP asymmetries of B — y K^p in both SM and the concerned new 
physics model under different scenarios. From the table, we firstly note that 5° — y Kq°p(co) have 
large asymmetries, and different scenarios have different signs but with large uncertainties. If we 
can calculate the annihilation accurately within some effective approach in future, this parameter 
could be used to distinguish the scenarios. Secondly, for B — » p°(co), the Z' could change 
the signs of the center values, and these two decays can be used in probing new physics effect. 

V. SUMMARY 

Motivated by recent measurements of decays B — >■ K^p and Kq§, we studied the branching 
fractions of these decays both in SM and in the family nonuniversal Z' model within the QCDF 
framework. Because it is not clear whether Kq is the lying state or the first orbital excited state, 
we calculate them under two different scenarios. For these decay modes with scalar meson, the 
weak annihilations play more important roles than that in B — y PP and PV decays, so that they will 
take large uncertainties. From this point, an effective way that could calculate the annihilations 
reliably is needed. Comparing with the experimental results, we found different channels favor 
different scenarios. Moreover, in order to account for the large isospin asymmetries in the data, 
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large weak annihilations are also required. Adding the contribution of the family nonuniversal Z ! 
boson, we note that both the branching fractions and their ratios are changed remarkably. However, 
we cannot identify the character of the scalar meson Kq, either. Furthermore, for most channels, 
the Z' contribution will be buried by large uncertainties, except for decays B — > Kq~ p (oo). 

In this work, we also calculated the CP asymmetries of these decays and found the CP asymme- 
tries of B — > Kq§ are almost zero. In different scenarios, the CP asymmetries of B — > Kq~p°(oo) 
have different signs, thus they can be used to classify the scalar Kq. If its character is identified, 
we accordingly could used these results to probe the new gauge boson Z', because it changes the 
signs of CP asymmetries. All above results could be tested in the running LHCb or the Super-b 
factories in future. 
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